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Abstract — We  review  the  salient  features  of  two  advanced 
nodes  of  an  8-Nb-layer  fully  planarized  process  developed 
recently  at  MIT  Lincoln  Laboratory  for  fabricating  Single  Flux 
Quantum  (SFQ)  digital  circuits  with  very  large  scale  integration 
(VLSI)  on  200-mm  wafers:  the  SFQ4ee  and  SFQ5ee  nodes,  where 
“ee”  denotes  the  process  is  tuned  for  energy  efficient  SFQ 
circuits.  The  former  has  eight  superconducting  layers  with 
0.5  pm  minimum  feature  size  and  a  2  fl/sq  Mo  layer  for  circuit 
resistors.  The  latter  has  nine  superconducting  layers:  eight  Nb 
wiring  layers  with  the  minimum  feature  size  of  350  nm  and  a  thin 
superconducting  MoNx  layer  ( Tc  ~  7.5  K)  with  high  kinetic 
inductance  (about  8  pH/sq)  for  forming  compact  inductors.  A 
nonsuperconducting  ( Tc  <  3  K)  MoNx  layer  with  lower  nitrogen 
content  is  used  for  6  f2/sq  planar  resistors  for  shunting  and 
biasing  of  Josephson  junctions  (JJ).  Another  resistive  layer  is 
added  to  form  interlayer,  sandwich-type  resistors  of  mf2  range 
for  releasing  unwanted  flux  quanta  from  superconducting  loops 
of  logic  cells.  Both  process  nodes  use  Au/Pt/Ti  contact 
metallization  for  chip  packaging.  The  technology  utilizes  one 
layer  of  Nb/A10x-Al/Nb  JJs  with  critical  current  density,  Jc  of 
100  pA/pm2  and  minimum  diameter  of  700  nm.  Circuit  patterns 
are  defined  by  248-nm  photolithography  and  high  density  plasma 
etching.  All  circuit  layers  are  fully  planarized  using  chemical 
mechanical  planarization  (CMP)  of  Si02  interlayer  dielectric. 
The  following  results  and  topic  are  presented  and  discussed:  the 
effect  of  surface  topography  under  the  JJs  on  the  their  properties 
and  repeatability,  7C  and  Jc  targeting,  effect  of  hydrogen  dissolved 
in  Nb,  MoNx  properties  for  the  resistor  layer  and  for  high  kinetic 
inductance  layer,  technology  of  mfl-range  resistors. 

Index  Terms — Josephson  junctions,  high  kinetic  inductance, 
MoNx,  Nb/A10x/Nb  Josephson  junctions,  ERSFQ,  RSFQ,  RQL, 
superconducting  integrated  circuits 


I.  Introduction 

SUPERCONDUCTING  DIGITAL  electronics  is  considered  for 
applications  in  high  performance-computing  due  to  a 
potential  for  much  higher  clock  rates  and  lower  energy 
dissipation  than  those  offered  by  the  current  CMOS 
technology.  In  order  to  become  competitive  with 
semiconductor  electronics,  superconducting  circuits  must 
reach  a  very  large  scale  of  integration  (VLSI)  that  would 
enable  circuit  functionalities  and  complexities  required  for 
computing.  At  present,  superconducting  digital  circuits  have 
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about  5  orders  of  magnitude  lower  integration  scale  than  the 
typical  CMOS  circuits.  E.g.,  the  largest  demonstrated  Single 
Flux  Quantum  (SFQ)  circuits  have  only  about  105  switching 
elements,  Josephson  junctions,  [l]-[3]  whereas  CMOS  circuits 
routinely  have  over  1010  transistors. 

There  are  several  reasons  for  this  gigantic  disparity: 
insufficient  funding  and  lack  of  profit-driven  investments  in 
superconductor  electronics  (SCE),  immaturity  of  the 
fabrication  processes  and  integrated  circuit  design  tools,  and 
physical  limitations  on  the  size  of  SFQ  logic  cells.  The  present 
work  is  devoted  to  development  of  advanced  fabrication 
processes  for  making  superconducting  VLSI  circuits  and,  thus, 
reducing  this  disparity. 

II.  Eight-Nb-Layer  Process  Node  SFQ4ee 

In  our  previous  works  [4], [5]  we  introduced  a  fabrication 
process  with  eight  Nb  layers  and  full  planarization  of  all  layers 
including  the  layer  of  Josephson  junctions.  This  process  node 
was  termed  SFQ4ee,  where  “ee”  denotes  that  the  process  is 
tuned  for  making  energy  efficient  circuits  for  the  IARPA 
Cryogenic  Computing  Complexity,  C3  Program  [6].  The  cross 
section  of  the  process  is  shown  in  Fig.  1;  the  target  parameters 
of  the  layers  as  they  appear  in  the  processing  and  minimum 
feature  sizes  (critical  dimensions)  are  given  in  Table  1. 
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a/c  is  the  minimum  critical  current  of  metal  lines  or  contact  holes  (vias)  of 
the  critical  dimensions  (CD)  in  mA,  Rs  is  the  sheet  resistance  of  resistor  layer 
R5  in  Q/sq. 

c  Josephson  critical  current  density  of  Nb/A10x-Al/Nb  trilayer  in  pA/pm2. 
d  Mixed  anodic  oxide  (AlNb)Ox  formed  by  anodization  of  Al/Nb  bilayer  of 
JJ  bottom  electrode. 


DF-4800  5.0kV  2.9mm  xl  0.0k  SE(U)  1 0/1 4/201 3 '  *  5.00um 


Fig.  1.  Scanning  electron  microscope  (SEM)  image  of  a  wafer  cross  section 
made  by  focus  ion  beam  (FIB).  The  wafer  was  fabricated  by  SFQ4ee  process. 
Labels  of  metal  layers  correspond  to  notations  in  Table  I.  Vias  corresponding 
to  the  physical  and  lithographical  layers  A#  and  I#  are  labeled  VO,. . ,,V6.  The 
vias  are  formed  by  etching  contact  holes  in  the  interlayer  dielectric  and 
covering  them  by  Nb  of  the  next  wiring  layer. 

A.  Fabrication  and  Measurement  Details 

All  metal  layers  (Nb,  Al,  Mo)  were  deposited  on  200  mm  Si 
wafers  by  magnetron  sputtering  using  a  multi-chamber  cluster 
tool  (Endura  from  Applied  Materials,  Inc.)  with  base  pressure 
of  10"8  Torr.  Si02  interlayer  dielectric  (ILD)  was  deposited  at 
150  °C,  using  a  Plasma  Enhanced  Chemical  Vapor  Deposition 
(PECVD)  system  Sequel  from  Novellus  (Lam  Research,  Inc.) 
Thickness  uniformity  of  the  deposited  oxide  was  a  =  2%, 
where  a  is  standard  deviation  (normalized  to  the  mean  value). 
Photolithography  was  done  using  a  Canon  FPA-3000  EX4 
stepper  with  248  nm  exposure  wavelength,  UV5  photoresist, 
and  AR3  bottom  antireflection  coating.  Etching  of  ah  metal 
and  dielectric  layers  was  done  in  a  Centura  etch  cluster 
(Applied  Materials,  Inc.),  using  Cl-based  chemistry  for  metals 
and  F-based  chemistry  for  dielectrics.  Planarization  of  the 
metal  layers  was  done  by  a  chemical  mechanical  planarization 
(CMP),  using  the  following  steps:  a)  deposition  of  a  -  2.5x 
times  thicker  Si02  over  the  patterned  metal  layer;  b)  polishing 
Si02  to  the  required  level,  using  a  CMP  tool  Mirra  from 
Applied  Materials,  Inc.  c)  measuring  the  remaining  dielectric 
thickness  in  49  points  on  the  wafer,  using  an  elipsometer;  d) 
redeposition  of  Si02,  if  needed,  to  achieve  the  target  ILD 
thickness  in  Table  1.  Etched  vias  10,  II,  etc.  are  filled  by  Nb  of 
the  following  metal  layer.  JJ  fabrication  was  described  in 
detail  in  [4]. 

For  the  process  characterization  and  electrical  testing  we 
used  a  Process  Control  Monitor  (PCM)  set  of  reticles  with 
sixteen  5  mm  x  5  mm  chips  containing  -  3600  structures  for 
testing  JJs,  wires,  vias,  inductors,  resistors,  dielectric  integrity, 
and  other  process  parameters.  For  purely  process  development 
fabrication  runs,  this  22  mm  x  22  mm  PCM  was  stepped  ah 
over  the  wafer  on  a  7  x  7  grid  using  labels  A,B,. .  .,G  for  the  x- 
coordinate  and  1,2,..., 7  for  the  y-coordinate.  For  runs  with 
SFQ  digital  circuits,  the  PCM  was  stepped  only  over  nine 
locations  (A4,  B2,  B6,  Dl,  D4,  D7,  F2,  F6,  G4)  whereas  a 
different  set  of  reticles,  containing  only  digital  circuits,  was 


used  for  the  other  40  locations. 

About  3300  test  structures  (TS)  per  each  location  were 
tested  at  room  temperature  by  using  a  semi-automated  wafer 
prober.  All  JJs  for  room- T  testing  were  designed  in  a  cross¬ 
bridge  Kelvin  resistor  (CBKR)  configuration  to  measure  the 
tunnel  barrier  resistance.  A  smaller  subset  of  TS,  including 
JJs,  wires,  vias,  resistors,  and  inductors,  were  tested  at  4.2  K 
after  dicing  the  wafer  and  wire  bonding  the  individual  chips.  A 
procedure  for  estimating  the  Josephson  critical  current  density 
Jc  and  JJ  parameter  spreads  from  room-  T  measurements  of  JJ 
conductance  was  described  in  [4].  Briefly,  a  junction 
conductance  G  at  room-T  is  directly  proportional  to  the  tunnel 
barrier  conductance  GN  at  4.2  K  and  to  /c  of  the  JJ  through  the 
Ambegaokar-Baratoff  relation  ( ICR N  product,  R^=  1/GN)  [7]. 
Accordingly,  the  relative  standard  deviation  of  conductances 
aG  (normalized  to  the  mean  conductance  value)  of  a  set  of 
nominally  identical  JJs  is  a  very  good  measure  of  the  relative 
standard  deviation  of  the  critical  current  aIc  in  the  set,  for 
details  see  [4]. 

Except  when  specifically  noted,  ah  JJs  were  designed  such 
that  Nb  leads  connecting  their  base  and  counter  electrodes  to 
100  pmx  100  pm  Au/Pt/Ti/Nb  contact  pads  were  interrupted 
by  50  Q  molybdenum  resistors,  25  pm  long  and  2  pm  wide. 
The  purpose  of  these  resistors  was  to  interrupt  a  diffusion  path 
from  the  JJs  to  contact  pads  through  Nb  wires  by  a  material 
with  much  lower  diffusion  coefficient  of  impurities  such  as 
hydrogen,  oxygen,  nitrogen,  etc.  These  impurities  may 
dissolve  in  Nb  during  various  steps  of  wafer  processing  and 
affect  junction  tunnel  conductance  and  Josephson  critical 
current  [8]-[12].  Molybdenum  resistors  work  as  a  diffusion 
barrier  preventing  out-diffusion  of  impurities  from  JJs  to 
contact  pads  and  thus  fixing  the  true  impurity  state  of  JJs  by 
the  end  of  wafer  processing  [10]-[12]. 

B.  Layer  Topography  Effects  on  Josephson  Junctions 

An  accurate  targeting  of  the  critical  currents  Ic  of  JJs  and 
minimization  of  their  variations  (standard  deviation  a/c)  are  of 
prime  importance  for  yielding  operational  SFQ  digital  circuits 
because  of  their  narrow  design  margins  with  respect  to 
deviations  of  /c  from  the  target  values.  One  of  the  possible 
sources  of  Ic  variations  is  surface  topography  created  by 
patterned  wiring  layers  under  JJs  as  a  result  of  planarization 
imperfections,  e.g.,  Si02  dishing  between  metal  lines,  etc.  For 
instance,  Fig.  2  -  Fig.  4  show  a  few  examples  of  JJs  placed 
over  the  edges  of  etched  wires  in  Nb  layers  below  the  JJs. 
Electric  properties  of  such  JJs  (tunnel  barrier  conductance  G  at 
room-  T)  have  been  measured  and  compared  with  a  control 
group  of  JJs  placed  over  the  flat  surface  formed  above  M4, 
M2,  and  M0  ground  planes  with  no  patterning  under  and  in 
proximity  to  the  JJs. 

As  a  ground  plane  we  refer  to  a  multiply-connected  metal 
layer  with  average  metal  density  of  >  80%  on  a  scale  of 
5  mm  x  5  mm  chips.  It  typically  serves  for  returning  ah  bias 
currents  in  circuits.  Since  electric  signals  coming  from  the 
bottom  layers  should  pass  though  the  ground  planes  in  order  to 
reach  JJs,  these  planes  have  multiple  etched  openings  (holes), 
moats,  and  short  wires  whose  edges  may  be  located  under  JJs, 
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e.g.,  as  shown  in  Fig.  2.  Wiring  layers  typical  have  25%  to  The  edges  of  patterns  in  Ml  and  M3  layers  were  about  50  pm 

65%  metal  density.  JJ  density  was  typically  from  4%  to  6%.  away  from  the  JJs  in  the  control  group. 


(a)  (b) 


Fig.  2.  SEM  images  of  the  FIB  cross  sections  through  a  JJ  and  Nb  layers 
below  it:  (a)  1  pm  circular  JJs  above  the  coincident  edge  of  wires  in  four  metal 
layers  MO,  Ml,  M2,  and  M4  below  the  junction;  (b)  a  zoom  in  showing  only 
the  JJs  and  the  edge  of  the  upper  layer  M4.  Layers  MO,  M2,  and  M4  were 
ground  planes  with  average  metal  density  ~  80%  on  the  chip  scale.  Please 
note,  FIB  samples  are  tilted  ~  55°,  so  the  vertical  scale  is  not  equal  to  the 
horizontal  scale. 

The  surface  of  the  14  dielectric  after  planarization  and  the 
surface  of  the  bottom  electrode  (M5)  of  the  JJ  look  very  flat  in 
Fig.  2,  despite  the  presence  of  the  coincident  edge  of  four 
wires  in  the  layers  below  right  in  the  middle  of  the  JJ.  Fig.  3 
gives  another  example  of  possible  topography  under  a  JJ. 
Here,  multiple  wires  of  the  minimum  allowed  width  (Table  1) 
are  placed  under  the  JJs.  Zoomed  image  Fig.  3(b)  does  not 
reveal  noticeable  changes  in  the  flatness  of  the  bottom 
electrode  of  the  JJ,  M5.  However,  the  tunnel  barrier  is  only  ~ 
1  nm  thick  and  thus  can  potentially  be  affected  by  surface 
imperfections  not  visible  at  this  magnification  or  in  these 
particular  cross-sections.  Therefore,  only  electric  testing  can 
reveal  if  the  presence  of  wire  edges  in  the  layers  under  the 
junction  layer  and  topography  they  create  affect  the  JJs. 


(a)  (b) 


Fig.  3.  (a)  FIB  cross  section  showing  a  1  pm  circular  JJ  above  topography 
created  by  multiple  0.5-pm-wide  M4  wires  stacked  on  0.5-pm-wide  M3  wires 
above  M2  and  MO  ground  planes;  (b)  a  zoom  in  showing  the  JJ  and  the  M4 
layer  under  it.  Globally,  M4  layer  was  a  ground  plane  with  the  metal  density 
about  80%.  The  M4  wires  shown  are  about  20  pm  long  in  the  perpendicular  to 
the  cross-section  direction  and  separated  from  the  ground  plane. 

Fig.  4  gives  an  example  when  the  planarization  does  not 
look  this  perfect.  In  this  case  a  few  0.7- pm- wide  parallel  wires 
in  layer  M3  were  placed  below  a  much  wider  wire  in  M4 
under  a  JJ.  Since  the  metal  density  of  the  M3  layer  near  the  JJ 
was  low,  there  was  a  noticeable  dishing  of  the  dielectric 
between  the  wires.  As  a  result,  M4  wires  sagged  a  little  bit 
between  them,  and  a  clear  step  in  M4  can  be  seen  projecting 
up  from  the  edge  of  the  M3  wire  right  under  the  JJ,  see 
Fig.  4(b). 

The  control  group  had  24  JJs  per  chip,  196  JJs  per  wafer, 
placed  above  M4,  M2,  and  M0  ground  planes,  as  shown  in 
Fig.  1  but  with  layers  Ml  and  M3  etched  away  under  the  JJs. 


. 

*41 


M4 
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(a)  (b) 

Fig.  4.  (a)  FIB  cross  section  showing  a  1  pm  circular  JJ  above  topography 
created  by  multiple  0.7-pm-wide  M3  wires  under  a  20-pm-wide  M4  wire  and 
above  M2  and  M0  ground  planes;  (b)  a  zoom  in  showing  the  JJ  and  the  M4 
layer  under  it.  Layer  M4  clearly  has  a  step  right  in  the  middle  of  the  JJ, 
resulting  from  sagging  in  (dishing  of  the  dielectric  between  two  wires  in  M3 
under  the  JJ.  It  appears  though  that  this  step  does  not  affect  the  planarity  of  the 
bottom  electrode  of  the  JJ,  but  only  electric  testing  can  reveal  if  this  JJ  differs 
in  any  way  from  a  JJ  placed  on  the  flat  surface  with  no  wires  below  it. 

Another  group  consisting  of  sixty  six  JJs  per  chip,  594  per 
wafer,  was  placed  over  the  topographies  shown  in  Fig.  2  - 
Fig.  4  and  a  few  other  combinations  of  layers  and  wires,  e.g., 
single  wire,  crossing  wires,  etc.  All  JJs  in  both  groups  were 
Manhattan- shaped  junctions  defined  on  a  5  nm  grid  with  an 
effective  diameter  of  1000  nm.  The  typical  aggregate 
distributions  of  conductance  of  the  JJs  in  both  groups  are 
shown  in  Fig.  5. 


Fig.  5.  Distribution  of  conductance  of  Josephson  junctions  on  a  wafer:  (a) 
junctions  above  ground  planes  M4,  M2  and  M0  with  no  patterning  under  and 
near  the  junctions,  (b)  junctions  above  various  patterns  in  M4  and/or  M3,  and 
M2  layers  described  in  the  text.  The  mean  values  <Gflt>  and  <Gtpgr>  as  well 
as  second  moments  of  the  distributions,  aflt  and  Otpgr,  are  shown.  Solid  lines 
show  fits  to  a  Gaussian  distribution.  The  distribution  of  conductance 


characterizes  the  distribution  of  the  critical  current.  Note  the  scale  difference 


between  (a)  and  (b). 
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From  the  point  of  view  of  SFQ  circuits,  the  important 
parameters  are  the  difference  in  the  mean  values  of  these 
distributions  <Gtpgr>  -  <Gfit>  and  the  difference  of  their 
second  moments  atpgr  -  opt,  where  subscripts  “tpgr”  and  “fit” 
refer  to  the  group  over  various  topographies  and  the  control 
group  above  the  nominally  flat  Si02  surface  above  the  ground 
plane,  respectively.  These  parameters  were  tracked  for  all 
wafers  made  by  the  process.  Hereafter  we  will  refer  to  the 
second  moments  of  conductance  and  Ic  distributions  as 
standard  deviations  without  implying  that  they  are  strictly 
Gaussian. 
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Fig.  6.  (a)  Relative  difference  of  the  mean  conductance  of  JJs  on  topographies 
<G!tpgr>  with  respect  to  the  junctions  on  the  flat  dielectric  above  M4  ground 
plane,  <Gflt>.  (b)  Difference  of  the  standard  deviations  of  these  two  groups  of 
JJs.  Parameters  of  aggregate  distribution  of  JJs  on  9  dies  per  wafer  are  shown 
by  squares  and  of  distributions  excluding  the  central  die,  D4,  are  shown  by 
stars. 


Fig.  6  shows  the  normalized  difference  of  the  mean  values 
(<Gtpgr>  -  <Gfit>)/<Gfit>  for  the  24  most  recently  made 
wafers.  It  can  be  seen  that  there  is  a  small  statistical  difference 
between  the  mean  values  of  conductance  of  JJs  over  various 
topographies  and  the  control  group,  from  about  0.5%  to  2.5%, 
and  a  somewhat  larger  difference  in  the  standard  deviation.  On 
average,  JJs  over  topographies  have  a  bit  higher  conductance 
and,  hence,  higher  critical  current  and  a  bit  wider  statistical 
distribution  than  JJs  above  the  ground  plane,  Fig.  6(b). 
However,  a  die-by-die  comparison  shows  that  the  main 
difference  is  confined  to  just  one  die,  the  central  die  D4,  on  all 


the  wafers.  Excluding  this  central  die  from  the  analysis  shows 
that  there  is  about  1%  increase  in  the  average  conductance  of 
JJs  over  topographies,  see  Fig.  6(a).  From  Fig.  6(b)  it  is  also 
clear  that  there  is  no  statistically  important  increase  in  the 
width  of  their  distribution,  given  that  the  estimated  uncertainty 
of  a  determination  using  JJ  resistance  measurements  on  a 
wafer  prober  is  about  ±0.3%. 

The  observed  increase  in  JJ  conductance  can  be  explained 
by  an  increase  in  JJ  average  area  resulting  from  purely 
geometric  factors  such  as  (macroscopic)  surface  roughness. 
For  instance,  assume  that  near  the  edge  of  metal  lines  shown 
in  Fig.  2(a)  and  Fig.  3(a)  a  step  h  in  the  dielectric  forms  as  a 
result  of  Si02  dishing  during  the  CMP,  as  shown 
schematically  in  Fig.  7(a).  For  small  steps  h  «  t,  where  t  is 
Nb  layer  thickness,  a  deposited  Nb/A10x-Al/Nb  trilayer  covers 
such  steps  likely  without  changes  in  the  tunnel  barrier 
properties.  The  step  causes  an  increase  of  the  junction  area 
with  respect  to  the  flat  JJ  of  about  SAtpgr  =  2 h(r2  -  a2)112,  where 
a  is  the  distance  between  the  step  and  the  center  of  the  circular 
junction  with  radius  r.  The  relative  increase  in  tunnel 
conductance  due  to  this  area  increase  is 
dG^g/Gfit  =  SAfpg/Af[X=  (2h/nr)(\  -  (a/rf)m. 


Fig.  7.  (a)  A  sketch  of  a  junction  placed  over  a  step  formed  in  the  interlayer 
dielectric  as  a  result  of  imperfect  planarization,  e.  g.,  dishing  of  Si02  near  the 
edges  of  wires  in  metal  layer  M3  and  M4  below  the  junction  base  electrode 
layer  M5.  A  bend  of  the  base  electrode  and  A10X-A1  barrier  with  height  h 
could  be  formed,  increasing  the  effective  area  of  the  junction,  (b)  Dishing  of 
interlayer  dielectric  between  two  metal  lines.  Due  to  the  surface  curvature,  a 
JJ  placed  between  these  lines  has  larger  area  than  the  one  on  the  flat  surface. 

The  observed  increase  in  the  mean  conductance  of  JJs  with 
r=  500  nm  is  ~  1%,  see  Fig.  6.  It  can  be  explained  by  the 
presence  of  steps  with  the  mean  height  of  about  8  nm.  They 
would  hardly  be  seen  in  SEM  images,  Fig.  3  and  Fig.  4.  The 
maximum  step  height  /zmax  we  observed  after  CMP  was  ~ 
35  nm,  which  could  result  in  the  maximum  relative  increase  in 
conductance  of  1000  nm  JJs  of  ~  4%.  We  note  that 
randomness  in  step  locations  and  step  heights,  i.e.,  the  surface 
roughness,  would  cause  random  variation  of  JJ  areas  and 
increase  the  width  of  JJ  conductance  distribution.  This 
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increase  of  the  normalized  aG  can  be  estimated  as  2ah/nr, 
where  ah  is  the  width  of  the  step  heights  distribution.  Then,  ah 
can  be  estimated  as  hmaJ6  ~  6  nm,  giving  aG  ~  4/r  about  0.4% 
for  1000  nm  JJs;  r  is  expressed  in  nm.  This  surface-roughness- 
caused  contribution  to  JJ  variability  is  of  the  same  magnitude 
as  the  contribution  of  the  so-called  mask  errors  estimated  in 
[4].  For  a  full  comparison  with  the  results  of  [4],  a  relationship 
between  the  actual  (on  wafer)  radius  of  JJs  r  and  the  designed 
radius  rd  should  be  used,  as  done  in  [4]. 

Another  example  is  Si02  dishing  between  two  metal  lines 
during  the  CMP,  forming  a  curved  surface  with  radius  R  and 
depth  /z,  see  Fig.  7(b).  A  JJ  formed  on  this  curved  surface  will 
have  larger  area  than  on  the  flat  surface.  Depending  on  the 
length  of  the  wires,  the  dishing  can  be  approximated  by  either 
a  cylindrical  or  a  spherical  surface  with  radius  R.  In  both  cases 
the  JJ  area  can  be  approximated  by  the  area  of  a  spherical 
sector  with  height  h,  At pgr  =  2nRk  =  n(r2  +  h2),  resulting  in  the 
increase  in  the  JJ  area  and  the  JJ  conductance  by  a  factor 
1  +(h/r)2.  For  1000  nm  JJs  and  h  =  35  nm,  the  increase  is 
about  0.5%. 

These  simple  estimates  demonstrate  how  macroscopic 
nonuniformities  (surface  roughness)  of  the  dielectric  surface 
under  the  JJs  may  affect  their  parameters,  even  without  taking 
into  account  possible  changes  in  the  microscopic  properties  of 
the  tunnel  barrier. 

JJs  located  near  the  very  center  of  200  mm  wafers  show,  on 
average,  larger  parameter  spreads  and  more  outliers  than  JJs 
located  more  than  ~  1cm  away  from  the  center.  Usually  these 
outliers  have  much  larger  conductance  than  the  mean, 
indicating  some  type  of  damage  to  the  tunnel  barrier  rather 
than  the  geometric  effects  considered  above.  E.g.,  such  an 
outlier  JJ  with  high  conductance  can  be  clearly  seen  in 
Fig.  5(b).  At  present  we  do  not  have  an  explanation  for  this 
behavior  and  its  cause. 

C.  Critical  Current  and  Jc  Targeting 

For  accurate  targeting  of  Ic  in  SFQ  circuits,  Josephson 
critical  current  density  Jc  is  a  convenient  design  parameter. 
However,  being  a  differential  characteristic,  its  accurate 
determination  requires  the  accurate  measurements  of  Ic  and  the 
accurate  knowledge  of  the  actual  size  of  JJs.  Also  the 
cryogenic  measurements  on  a  200  mm  wafer  scale  are 
somewhat  involved.  The  full  procedure  was  described  in  [4]. 
For  the  purpose  of  this  paper,  we  looked  at  the  targeting  of  the 
critical  current  (tunnel  conductance)  of  JJs  with  design 
diameter  of  1000  nm  because  it  is  perhaps  the  most  frequently 
used  size  of  JJs. 

Fig.  8  shows  the  wafer-averaged  conductance  of  1000  nm 
JJs  on  64  wafers  fabricated  from  the  beginning  of  2014  until 
August  2015.  A  set  of  216  JJs  per  wafer  at  the  nine  PCM 
locations  described  before  was  used.  The  second  moment  of 
the  distributions  (standard  deviation)  aG  of  conductance  in  the 
group  is  shown  as  an  error  bar.  For  JJs  of  this  diameter,  tunnel 
barrier  conductance  of  50  mS  corresponds  approximately  to 
the  target  Jc  of  100  pA/pm2. 

A  ±10%  margin  band  is  shown  by  the  dashed  lines,  which  is 
assumed  to  be  the  acceptable  margin  for  SFQ  circuits.  It  is 
easy  to  count  that  33  out  of  the  64  studied  wafers  were  within 


this  margin,  giving  a  Jc  targeting  wafer  yield  Y  of  52%.  Since 
the  control  of  the  JJs  area  and  its  run-to-run  repeatability  in 
our  process  is  very  tight,  the  changes  in  JJ  average 
conductance  in  Fig.  8  are  fully  attributable  to  changes  in  the 
tunnel  barrier  transparency  and  Jc  from  run  to  run.  We 
continue  investigating  the  factors  affecting  Jc  targeting  in  our 
process,  beyond  the  oxidation  time  and  pressure,  and  how  this 
targeting  can  be  significantly  improved. 


Fig.  8.  Tracking  of  the  wafer-averaged  conductance  of  junctions  at  room 
temperature.  JJs  were  designed  as  manhattans  on  a  5-nm  grid  and  were 
circularly  shaped  with  the  effective  diameter  of  1000  nm.  Red  circles 
correspond  to  the  JJs  with  molybdenum  resistors  in  Nb  wires  leading  to  the 
contact  pads  and  black  squares  to  the  JJs  without  the  resistors.  Mo  resistors 
work  as  diffusion  barriers  and  do  not  allow  hydrogen  dissolved  in  the  base 
and  counter  electrodes  of  JJs  and  in  the  top  Nb  wire  to  escape  into  the  contact 
pads.  The  changes  observed  reflect  the  run-to-run  variation  of  the  tunnel 
barrier  transparency  and  Jc. 

D.  Hydrogen  Effect  on  Critical  Current  and  Jc  Targeting 

It  is  well  known  that  many  impurities,  especially  hydrogen, 
can  dissolve  in  Nb  during  wafer  processing  and  reduce  Jc  and 
conductance  of  JJs  due  to  their  effect  on  the  tunnel  barrier 
height  [10], [12].  A  convenient  method  of  monitoring  the 
degree  of  hydrogen  contamination  of  JJs  is  based  on 
comparing  two  groups  of  JJs:  JJs  directly  connected  to 
Au/Pt/Ti/Nb  contact  pads  by  Nb  wires  and  JJ  connected 
through  molybdenum  resistors,  which  are  a  diffusion  barrier 
for  hydrogen.  Hydrogen  can  escape  from  the  JJs  to  the  contact 
pads  in  the  first  group  but  stays  inside  the  JJs  in  the  second 
group.  This  leads  to  a  higher  conductance  and  higher  Jc  of  JJs 
in  the  first  group  than  in  the  second  [12].  The  diffusion  path 
from  the  JJs  to  the  contact  pads  is  short,  no  more  than  50  pm. 
So  hydrogen  out-diffusion  happens  during  a  very  short  time 
period  after  the  contact  pad  deposition  [12]-[14]. 

Fig.  8  shows  the  wafer-averaged  conductance  for  these  two 
groups  of  JJs  for  64  wafers  studied.  We  use  notations  <Gwr> 
and  <Gr>  for  the  mean  conductance  of  the  group  without  the 
resistors  and  the  group  with  resistors,  respectively.  Their 
normalized  difference  (<Gwr>-<Gr>)/<Gr>  gives  the  relative 
difference  in  the  mean  values  of  Ic  and  Jc  in  these  two  groups 
and  is  shown  in  Fig.  9.  Knowing  this  difference  may  be 
important  for  designing  circuits  which  have  junctions  directly 
connected  to  the  contact  pads  by  their  counter  electrodes  or 
other  types  of  connections  that  may  create  a  difference  in  the 
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hydrogen  content  and  hence  create  bi-modal  distributions  of 
nominally  identical  JJs  within  the  circuit. 


0  10  20  30  40  50  60  70 


Wafer  # 


Fig.  9.  The  relative  difference  of  JJ  conductance  of  junctions  without 
hydrogen  dissolved  in  Nb  (without  molybdenum  diffusion  barrier),  Gm,  and 
with  hydrogen  dissolved  in  Nb  during  the  processing  (and  preserved  by  the 
presence  of  molybdenum  resistors),  Gw.  This  quantity  gives  also  the  relative 
difference  of  the  critical  currents  of  JJs  in  these  two  groups:  without  dissolved 
hydrogen  in  Nb  and  with  dissolved  hydrogen.  JJ  conductance  data  are  wafer- 
averaged.  The  same  set  of  wafers  as  in  Fig.  8  was  used. 


e)  an  additional  resistive  layer  is  added  between  Nb  layers 
M4  and  M5  in  order  to  enable  interlayer,  sandwich-type 
resistors  with  resistance  values  in  the  mQ  range  for 
minimizing  magnetic  flux  trapping  and  releasing  unwanted 
flux  from  logic  cells. 


Dn<'i  o+Ar  DC 


Layer 


Fig.  10.  Cross  section  of  a  wafer  fabricated  by  the  SFQ5ee  process.  The  labels 
of  metal  layers  and  vias  are  the  same  as  in  Table  I.  New  features  of  the 
SFQ5ee  are  shown:  a  high  kinetic  inductance  layer  under  M0  and  a  layer  of 
mD-range  resistors  between  M4  and  M5  layers. 


It  can  be  seen  from  Fig.  9  that  the  level  of  hydrogen 
contamination  varies  from  wafer  to  wafer.  We  have  many 
indications  that  the  source  of  hydrogen  contamination  is 
PECVD  interlayer  dielectric  which,  due  to  the  low  deposition 
temperature,  has  high  concentration  of  Si-OH  and  Si-H  bonds 
and  can  also  absorb  H20  molecules  during  the  processing. 
Some  of  the  loosely  bonded  hydrogen  can  diffuse  into  Nb 
during  processing  steps  taking  place  at  elevated  temperatures. 
However,  this  needs  further  investigating  as  well  as  finding 
the  methods  of  eliminating  this  contamination.  We  note  here 
that  hydrogen  diffusion  barriers  proposed  and  implemented  in 
[11]  do  not  eliminate  hydrogen  contamination  of  Nb  but 
simply  fix  the  hydrogen  content  by  preventing  hydrogen 
migration  between  the  circuit  layers.  As  can  be  seen,  the 
difference  between  the  JJs  with  hydrogen  and  without 
hydrogen  can  reach  ~  10%.  Therefore,  we  continue 

investigating  all  possible  ways  of  minimizing  hydrogen 
contamination  of  Nb  in  our  process. 

III.  SFQ5ee  Process  Node 

In  comparison  with  the  SFQ4ee  features  listed  in  Table  1,  a 
more  advanced  SFQ5ee  node  offers  the  following 
improvements: 

a)  the  minimum  linewidth  and  spacing  for  all  metal  layers, 
except  M0  and  R5,  is  reduced  to  0.35  pm  and  0.5  pm, 
respectively; 

b)  the  minimum  size  of  etched  vias  and  their  metal  surround 
is  reduced  to  0.5  pm  and  0.35  pm,  respectively; 

c)  the  sheet  resistance  of  the  resistor  layer  is  increased  to 
6  Q/sq  by  utilizing  a  nonsuperconducting  MoNx  film,  offering 
a  choice  of  either  2  Q/sq  (Mo)  or  6  f2/sq  (MoNx)  planar 
resistors  for  JJ  shunting  and  biasing; 

d)  an  additional  thin  superconducting  layer  with  high 
kinetic  inductance  is  added  below  the  first  Nb  layer  M0  in 
order  to  enable  compact  bias  inductors; 


A  cross-section  showing  these  new  features  of  the  SFQ5ee 
node  is  given  in  Fig.  10.  Below  we  will  discuss  new  features 
c),  d),  and  e)  and  the  motivation  for  having  them  in  some 
detail,  leaving  the  details  of  the  fabrication  processes  involved 
and  the  detailed  results  of  parameter  characterization  for  a 
separate  publication. 

A.  High-Shee-  Resistance  Layer 

Shunt  resistors  in  the  SFQ4ee  process  occupy  a 
considerably  larger  area  than  JJs,  restricting  circuit  integration 
scale.  It  may  seem  easy  to  decrease  the  area  of  resistors  by 
increasing  the  sheet  resistance  of  the  resistor  layer.  Indeed,  a 
thin-film  resistor  with  R  =  ( l/w)Rs  occupies  area 
Ar  =  w2R/Rs  +  2 Ac  where  Rs  is  the  sheet  resistance,  w  and  /  are 
the  resistor  width  and  length  between  two  superconducting 
contacts  through  the  dielectric  (vias),  and  2^4C  is  the  area  of 
these  contacts,  see  the  top  of  Fig.  11.  So,  increasing  Rs  should 
decrease  AR. 

The  expression  above  assumes  that  the  length  /  can  be  made 
arbitrarily  small  in  order  to  get  small  values  of  R.  However, 
this  is  not  the  case.  The  superconducting  vias  to  the  resistor 
cannot  be  placed  arbitrarily  close  to  each  other.  They  should 
be  separated  by  the  minimum  distance  /min  equal  to  the 
minimum  spacing,  s  between  Nb  wires  plus  twice  the 
minimum  required  surround,  sr  of  the  via  by  the  wire.  In  the 
SFQ4ee  process,  sr  =  0.5  pm  and  s=  0.7  pm,  making 
/min  =1.7  pm.  This  sets  the  minimum  value  of  the  resistor  that 
can  be  made  using  the  minimum  allowed  width  in  the  process, 
wmin,  and  hence,  with  the  minimum  possible  area  as  R  = 
(/min/jvmin)^s  =  6.8  Q  in  the  SFQ4ee  node.  For  getting  the 
smaller  values,  the  resistor  length  cannot  be  reduced  below 
/m1T1.  Therefore,  the  width  needs  to  be  increased  with  respect  to 
the  minimum  allowed.  In  this  regime,  the  resistor  area  is  given 
by  Ar  =  (lmin)2Rs/R  +  2 Ac  and  increases  with  increasing  Rs. 

Fig.  11  shows  dependences  of  the  resistor  active  area 
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(without  the  area  of  vias)  as  a  function  the  resistor  value  for 
Rs  =  2  Q/sq  and  6  Q/sq.  The  minimum  achievable  area 
(/minWmin)  does  not  depend  on  Rs  and  is  reached  at 

R  (4nin/WminV^s« 


u 
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Fig.  11.  Shunt  resistor  area  (without  the  area  of  Nb  contacts)  for  different 
values  of  the  resistor:  SFQ4ee  process  with  Rs  =  2  Q/sq,  wmin=0.5  pm,  and 
/min  =1.7  pm  (black  solid  curve),  and  SFQ5ee  with  Rs  =  6  Q/sq,  wmjn=  0.5  pm, 
and  /min=  1.2  pm  (top,  blue  solid  curve).  The  regime  shown  by  the  dashed 
lines  cannot  be  realized  because  the  length  /  of  the  resistor  cannot  be 
decreased  below  /min.  The  doted  curve  corresponds  to  a  future  SFQ6ee  process 
with  Rs  =  6  Q/sq,  wmin  =  s  =  0.35  pm  and  =  0.95  pm.  The  top  scale  shows 
the  area  of  the  critically  damped  junction  with  shunt  resistor  value  R.  A  sketch 
of  the  resistor  cross  section  is  shown  on  top. 

For  the  Jc  =  100  pA/pm2  processes,  the  critical  damping  of 
JJs  corresponds  to  the  characteristic  voltage  Vc  =  IcRn  of 
-700  pV,  where  Rn  is  a  parallel  combination  of  the  shunt 
resistor  R  and  junction  subgap  resistance  7?sub.  The  latter  is 
about  lOx  times  larger  than  the  junction  normal-state 
resistance  R N  and  can  be  neglected  for  the  estimates  here. 
Then,  the  shunt  resistor  (in  ohms)  is  given  approximately  by 
R  =  1/Ah  where  Aj  is  the  junction  area  in  pm2.  Hence,  all  JJs 
with  areas  Aj  >  7/(3 Rs)  have  shunts  limited  by  the  minimum 
spacing  and  therefore  increasing  the  sheet  resistance  of  the 
shunt  material  without  reducing  /min  will  only  increase  the 
resistor  area.  For  the  SFQ4ee  process  this  corresponds  to  Aj 
>1.17  pm2  ( Ic  >117  pA),  covering  nearly  the  whole  range  of 
JJs  used  in  RSFQ  [15]  and  ERSFQ  [16]  circuits. 

If  we  increase  Rs  to  6  Q/sq,  keep  the  same  minimum  wire 
spacing  s  =  0.5  pm  and  reduce  the  surround  to  sr  =  0.35  pm  as 
in  the  SFQ5ee  process,  all  resistors  with  values  below  8.31  Q 
will  still  increase  in  area.  This  will  affect  resistive  shunts  of  all 
JJs  with  A j  >  0.84  pm2,  /c  >  84  pA,  i.e.,  all  junctions  used  in 


RSFQ  and  ERSFQ  circuits.  Other  types  of  logic,  which  use 
smaller  JJs  and  larger  shunts,  e.g.,  RQF  [17]  may  still  benefit 
from  the  increase  in  Rs.  If  in  the  next  process  node  we  also 
reduce  the  minimum  spacing  to  0.35  pm,  the  shunt  area 
savings  will  take  place  for  the  entire  range  of  resistor  values 
and  JJs  used  in  SFQ  circuits,  as  shown  by  a  dotted  line  in 
Fig.  11. 

Despite  the  absence  of  resistor  area  savings  in  SFQ5ee  node 
for  resistors  with  R<  8.31  Q,  we  have  developed  MoNx 
resistors  with  the  target  Rs  =  6  Q/sq  for  the  use  in  the  future 
process  nodes  with  smaller  resistor  linewidth  and  wire  spacing 
wmin  =  s  =  0.35  pm  and  0.25  pm,  and  also  for  the  nodes  with 
higher  Jc  such  as  200  pA/pm2  and  500  pA/pm2.  Higher  Jc 
junctions  require  larger  shunt  resistors  and  increasing  Rs  may 
be  beneficial. 

Amorphous  MoNx  films  with  a  broad  range  of  sheet 
resistances  and  superconducting  critical  temperatures  were 
obtained  by  changing  N2/Ar  ratio  during  Mo  sputtering  [18]- 
[20].  Fig.  12  shows  a  dependence  of  Tc  on  Rs  for  amorphous 
MoNx  films  with  40  nm  thickness,  deposited  on  thermally 
oxidized  Si  wafers  at  room-T.  For  the  Rs  range  of  interest,  the 
resistance  ratio  RRR  =  R300/R4.2  is  from  1.2  to  1.3,  so  we  target 
the  films  with  about  7.5  Q/sq  at  300  K  for  the  shunt  resistor 
layer  in  the  SFQ5ee  process. 
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Fig.  12.  Critical  temperature  of  amorphous  MoNx  films  deposited  on  Si02/Si 
wafers  at  room  temperature  by  reactive  sputtering  of  Mo  in  Ar/N2  mixture. 
Films  with  Tc  below  about  3.5  K  can  be  used  for  the  high  sheet  resistance 
resistor  layer.  High-Tc  (high  sheet  resistance)  region  can  be  used  for  a  layer  of 
inductors  with  high  kinetic  inductance  Lk  in  the  range  from  about  4  pH/sq  to 
about  16  pH/sq  at  40  nm  film  thickness.  Dashed  line  is  to  guide  the  eye  only. 


Sheet  resistance  uniformity  of  unpattemed  MoNx  films  of 
better  than  2%  on  200-mm  wafers  was  achieved  by  rotating 
the  wafer  during  the  film  deposition.  Four-probe 
measurements  were  done  in  49  points.  On  the  fully  processed 
SFQ5ee  wafers,  the  patterned  resistor  uniformity  (la  standard 
deviation)  was  measured  using  9  PCM  locations  and  found  to 
be  2.3%. 

B.  High-Kinetic-  Inductance  Layer 

It  can  be  easily  shown  that  inductors  occupy  most  of  the 
space  in  SFQ  circuits  of  any  type  RSFQ,  RQF,  AQFP,  etc. 
Making  energy  efficient  ERSFQ  circuits  requires  particularly 
large,  -  1 00  pH,  bias  inductors  instead  of  bias  resistors  in 
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order  to  eliminate  the  static  power  dissipation  [16].  A  100  pH 
inductor  in  the  SFQ4ee  process  node  occupies  more  than 
100  pm2  area,  which  severely  limits  the  integration  scale  of 
ERSFQ  circuits.  In  order  to  alleviate  restrictions  related  to  the 
size  of  inductors  and  enable  increasing  the  integration  scale  of 
ERSFQ  circuits,  we  suggested  and  implemented  in  SFQ5ee 
node  an  additional  (the  9th)  superconducting  layer  with  high 
kinetic  inductance.  It  is  well  known  that  kinetic  inductance  of 
thin  superconducting  films  Lk  =  fi^/t  is  much  larger  that  their 
geometric  inductance  if  X  »  t  and  X2/t»w ,  where  X  is 
magnetic  field  penetration  depth,  t  and  w  are  the  film  thickness 
and  width,  respectively. 

Many  high-kinetic-inductance  materials  such  as  NbNx, 
TiNx,  and  TiSi2  have  been  extensively  studied  for  applications 
in  detectors  of  photons  and  particles  [21].  Usually  such 
detectors  require  only  one  superconducting  layer  and  operate 
below  the  FHe  temperature,  so  very  thin  films  with  t  ~  5  nm 
can  be  used  giving  very  large  values  of  Lk.  In  our  case  we 
needed  to  integrate  a  thin  superconducting  film  into  a 
multilayer  process  and  prevent  a  significant  degradation  of  its 
critical  temperature  Tc  in  comparison  with  Tc  of  Nb.  Also  we 
need  to  keep  the  deposition  temperature  relatively  low  and 
would  like  to  minimize  the  number  of  different  materials  used 
in  the  process  because  each  additional  material  requires  an 
additional  vacuum  chamber  and  a  sputtering  gun. 

In  SFQ4ee  process  we  use  pure  molybdenum  thin  films  for 
the  2  fl/sq  resistor  layer.  In  SFQ5ee  process  we  add  nitrogen 
impurities  to  molybdenum  to  form  an  amorphous  MoNx  with 
low  nitrogen  content  [20]  in  order  to  increase  the  sheet 
resistance  to  6  fl/sq.  Therefore,  we  decided  to  investigate  and 
implement  thin  superconducting  MoNx  films  in  the  high-sheet- 
resistance  region  of  Fig.  12  (with  high  nitrogen  content)  as 
high  kinetic  inductors  [22].  The  advantage  of  this  approach  is 
that  the  same  vacuum  chamber  and  Mo  target  can  be  used  for 
reactive  sputtering  of  both  the  high- sheet-resistance  layer  for 
resistors  and  the  high-kinetic-inductance  layer  for  inductors 
simply  by  adjusting  the  N2/Ar  ratio  in  the  mixture. 

The  high-kinetic-inductance  layer  was  placed  as  the  first 
layer  in  the  process  as  shown  in  Fig.  10.  We  used  a  40  nm 
thick  MoNx  film.  After  patterning  it  was  covered  by  a  60  nm 
layer  of  Si02  to  isolate  it  from  the  first  Nb  layer  M0. 
Superconducting  vias  were  made  by  etching  contact  holes  in 
Si02.  Due  to  its  low  thickness,  the  layer  of  kinetic  inductors 
was  not  planarized  by  CMP. 

Sheet  inductance  of  the  unpattemed  films  was  measured 
using  a  dielectric  resonator  technique  [23], [24].  On  the  fully 
processed  wafers,  kinetic  inductors  of  various  widths  from 
0.5  pm  to  2  pm  and  various  lengths  were  included  in  the  arms 
of  the  test  SQUIDs  described  in  [5].  The  inductance  was 
determined  from  the  period  of  the  SQUID  modulation.  Sheet 
inductance  values  in  the  range  from  about  4  pH/sq  to  1 6  pH/sq 
for  40 -nm- thick  films  were  obtained  by  varying  the  MoNx 
deposition  parameters.  Higher-kinetic-inductance  films  have 
lower  superconducting  critical  currents.  The  target  value  for 
the  kinetic  inductance  in  the  SFQ5ee  process  node  was  set  as 
8  pH/sq  in  order  to  satisfy  the  critical  current  requirements  for 
bias  inductors  in  ERSFQ  circuits,  Ic>  0.5  mA  per  inductor. 


The  details  of  MoNx  high  kinetic  inductor  fabrication  and 
characterization  will  be  given  elsewhere. 

C.  Interlayer,  Sandwich-Type  mQ-Range  Resistors 

One  of  the  main  reliability  problems  of  SFQ  circuits  is  the 
sensitivity  of  their  operation  margins  to  magnetic  flux  trapping 
inside  and  near  the  logic  cells.  The  flux  is  caused  by  the 
residual  magnetic  field  and  the  field  produced  by  bias 
currents.  Single  flux  quanta  are  also  produced  as  a  result  of 
cells  normal  operation.  If  trapped  in  superconducting  loops, 
some  types  of  cells  may  not  operate  correctly  as  a  result.  In 
order  to  minimize  both  types  of  the  flux  trapping,  it  would  be 
helpful  to  break  some  of  the  superconducting  loops  by 
resistors  allowing  magnetic  flux  to  decay.  The  typical  decay 
time  t  is  L/R ,  where  L  and  R  are  the  loop  inductance  and  the 
breaking  resistance.  Most  of  the  inductors  in  SFQ  circuits  are 
on  the  order  of  10  pH.  The  decay  time  should  be  much  longer 
than  the  circuit  clock  period  l/fc\  to  allow  SFQ  information 
processing.  On  the  other  hand  it  should  be  short  enough  to 
prevent  memory  effects,  i.e.,  to  reset  the  cell  to  a  zero  flux 
state  before  the  next  logic  operation.  At  1  GHz  clock 
frequency,  these  requirements  result  in  0  «  R  «  0.1  O. 

Due  to  the  limitations  described  in  III.  A,  it  is  not  possible  to 
make  planar  resistors  with  such  small  resistance  values  and 
reasonably  small  area  by  using  conventional  thin  films  having 
~  1  fl/sq  sheet  resistances.  However,  if  a  resistive  film  is 
sandwiched  between  two  superconducting  layers,  its  resistance 
in  the  direction  perpendicular  to  the  film  plane  can  be  utilized, 
giving  R  =  Rst2/A,  where  t  is  the  film  thickness  and  A  is  the 
contact  area.  For  any  t  <  0.1  pm  and  Rs  <  10  Q/sq,  RA  is  less 
than  0.1  Qpm2  and  the  required  compact  resistors  can  be 
realized.  The  main  trick  is  to  prevent  the  formation  in  this 
manner  of  a  sandwich-type  SNS  Josephson  junction  with 
critical  current  larger  than  O0/T. 

In  this  work,  a  resistive  film  was  sandwiched  between  Nb 
layers  M4  and  M5  as  shown  in  Fig.  10.  In  order  to  achieve 
this,  the  vias  between  these  layers  were  made  as  Nb  pillars 
(studs),  as  can  be  seen  in  Fig.  10,  by  using  a  double-etch  and 
planarization  (DEAP)  technology  of  stud-vias  described  in 
[25], [26].  After  the  CMP  to  the  level  of  Nb  pillars,  the 
resistive  layer  was  deposited  over  and  patterned  to  leave  it 
only  on  the  studs  were  a  resistive  contact  is  required.  After 
that  the  next  Nb  layer  M5  was  deposited,  completing  the 
sandwich.  In  this  structure,  the  contact  area  A  is  determined  by 
the  area  of  the  Nb  pillar  (stud- via).  More  details  on  the 
materials  used,  processing,  and  characterization  of  vertical, 
sandwich-type  resistors  will  be  given  elsewhere  due  to  a  lack 
of  space  here. 

Acknowledgment 

The  authors  would  like  to  thank  Vasili  K.  Semenov,  Alex  F. 
Kirichenko,  Marc  A.  Manheimer,  and  D.  Scott  Holmes  for 
useful  discussions,  Scott  Zarr  and  Corey  Stull  for  help  in  the 
SFQ5ee  process  development. 


2A-E-01.1 


9 


References 

[1]  S.  Nagasawa,  K.  Hinode,  T.  Satoh,  M.  Hidaka,  H.  Akaike,  A.  Fujimaki, 
N.  Yoshikawa,  K.  Takagi,  N.  Takagi,  “Nb  9-layer  fabrication  process  for 
superconducting  large-scale  SFQ  circuits  and  its  process  evaluation,” 
IEICE  Trans.  Electron.,  vol.  E97-C,  No.  3,  pp.  132-140,  March  2014. 

[2]  M.  Hidaka,  S.  Nagasawa,  K.  Hinode,  and  T.  Satoh,  “Device  yield  in  Nb- 
nine-layer  circuit  fabrication  process,”  IEEE  Trans.  Appl.  Supercond., 
vol.  23,  No.  3,  1100906,  June  2013. 

[3]  V.K.  Semenov,  Yu.A.  Polyakov,  and  S.K.  Tolpygo,  “New  AC-powered 
SFQ  digital  circuits,”  IEEE  Trans.  Appl.  Supercond.,  vol.  25,  No.  3, 
1301507,  June  2015. 

[4]  S.K.  Tolpygo,  V.  Bolkhovsky,  T.J.  Weir,  L.M.  Johnson,  M.A.  Gouker, 
and  W.D.  Oliver,  “Fabrication  process  and  properties  of  fully-planarized 
deep-submicron  Nb/Al-A10x/Nb  Josephson  junctions  for  VFSI  circuits,” 
IEEE  Trans.  Appl.  Supercond.,  vol.  25,  No.  3,  1 101312,  June  2015. 

[5]  S.K.  Tolpygo,  V.  Bolkhovsky,  T.J.  Weir,  C.J.  Galbraith,  F.M.  Johnson, 
M.A.  Gouker,  and  V.K.  Semenov,  “Inductance  of  circuit  structures  for 
MIT  FF  superconductor  electronics  fabrication  process  with  8  niobium 
layers,”  IEEE  Trans.  Appl.  Supercond.,  vol.  25,  No.  3,  1100905,  June 
2015. 

[6]  M.A.  Manheimer,  “Cryogenic  computing  complexity  program:  Phase  1 
introduction,”  IEEE  Trans.  Applied.  Supercond.,  vol.  25,  No.  3, 
1301704,  June  2015. 

[7]  V.  Ambegaokar  and  A.  Baratoff,  “Tunneling  between  superconductors,” 
Phys.  Rev.  Lett.,  vol.  10,  no.  11,  pp.  486^189,  June  1963. 

[8]  K.  Hinode,  T.  Satoh,  S.  Nagasawa,  and  M.  Hidaka,  “Hydrogen- 
inclusion-induced  variation  of  critical  current  in  Nb/AlOx/Nb  Josephson 
junctions,”  J.  Appl.  Phys.,  vol.  104,  no.  2,  July  2008,  Art.  ID  023909. 

[9]  S.K.  Tolpygo  and  D.  Amparo,  “Fabrication-process-induced  variations 
of  Nb/Al/AlOx/Nb  Josephson  junctions  in  superconductor  integrated 
circuits,”  Supercond.  Sci.  Technol.,  vol.  23,  No.  2,  034024,  Feb.  2010. 

[10]  K.  Hinode,  T.  Satoh,  S.  Nagasawa,  and  M.  Hidaka,  “Origin  of  hydrogen- 
inclusion-induced  critical  current  deviation  in  Nb/AlOx/Nb  Josephson 
junctions,”/.  Appl.  Phys.,  vol.  107,  No.  7,  073906,  Apr.  2010. 

[11]  S.K.  Tolpygo,  D.  Amparo,  R.T.  Hunt,  J.A.  Vivalda,  and  D.T.  Yohannes, 
“Diffusion  stop-layers  for  superconducting  integrated  circuits  and  qubits 
with  Nb-based  Josephson  junctions,”  IEEE  Trans.  Appl.  Supercond., 
vol.  21,  No.  3,  pp.  1 19-125,  June  2011. 

[12]  D.  Amparo  and  S.K.  Tolpygo,  “Investigation  of  the  role  of  H  in 
fabrication-process-induced  variations  of  Nb/Al/AlOx/Nb  Josephson 
junctions,”  IEEE  Trans.  Appl.  Supercond.,  vol.  21,  No.  3,  pp.  126-130, 
June  2011. 

[13]  J.  Engelhard,  “The  diffusion  of  H  and  D  in  Nb  and  Ta  at  low 
temperatures,”  J.  Phys.  F,  vol.  9,  no.  11,  pp.  2217-2229,  1979. 

[14]  J.H.  Scofield  and  W.W.  Webb,  “Resistance  fluctuations  due  to  hydrogen 
diffusion  in  niobium  thin  films,”  Phys.  Rev.  Lett.,  vol.  54,  pp. 353-356, 
Jan.  1985. 

[15]  K.K.  Fikharev  and  V.  K.  Semenov,  “RSFQ  logic/memory  family:  A  new 
Josephson-junction  technology  for  sub-terahertz-clock-frequency  digital 
systems,”  IEEE  Trans.  Appl.  Supercond.,  vol.  1,  no.  1,  pp.  3-28,  Mar. 
1991. 

[16]  D.F.  Kirichenko,  S.  Sarwana,  and  A.F.  Kirichenko,  “Zero  static  power 
dissipation  biasing  of  RSFQ  circuits,”  IEEE  Trans.  Appl.  Supercond., 
vol.  21,  no.  3,  pp.  776-779,  June  2011. 

[17]  Q.P.  Herr,  A.Y.  Herr,  O.T.  Oberg,  and  A.G.  Ioannidis,  “Ultra-low-power 
superconducting  logic,”  J.  Appl.  Phys.,  vol.  109,  2011,  Art.  ID.  103903. 

[18]  H.  Yamamoto,  T.  Miki,  M.  Tanaka,  “Preparation  of  superconducting 
MoNx  films  by  reactive  sputtering,”  in  Advances  in  Cryogenic 
Engineering  Materials,  vol.  32,  ch.  80,  R.P  Reed  and  A.F  Clark,  Eds. 
Springer,  1986,  pp.  671-678. 

[19]  Y.H.  Shi,  B.R.  Zhao,  Y.Y.  Zhao,  F.  Fi,  and  J.R.  Fiu,  “Superconducting 
and  normal-state  properties  of  MoNx  thin  films,”  Phys.  Rev.  B,  vol.  7, 
pp.  4488-4491,  Sep.  1988. 

[20]  Yohannes,  R.T.  Hunt,  J.A.  Vivalda,  D.  Donnelly,  D.  Amparo,  and 
A.F.  Kirichenko,  “20  kA/cm2  process  development  for  superconducting 
integrated  circuits  with  80  GHz  clock  frequency,”  IEEE  Trans.  Appl. 
Supercond.,  vol.  17,  No.  2,  pp.  946-951,  June  2007. 

[21]  J.  Zmuidzinar,  “Superconducting  microresonators:  physics  and 
applications,” /wzw.  Rev.  Cond.  Matt.  Phys.,  vol.  3,  pp. 169-214,  2012. 

[22]  A.  Patel,  A.-D.  Brown,  W.-T.  Hsieh  et  al.,  “Fabrication  of  MKIDS  for 
the  MicroSpec  spectrometer,”  IEEE  Trans.  Appl.  Supercond.,  vol.  23, 
no.  3,  June  2013,  Art.  ID  2400404. 

[23]  Z.Y.  Shen,  C.  Wlker,  P.  Pang,  W.F.  Holstein,  D.  Face,  and  D.J.  Kountz, 
“High  Tc  superconductor-sapphire  microwave  resonators  with  extremely 


high  Q-values  up  to  90  K,”  IEEE  Trans.  Microwave  Theor.  Tech., 
vol.  40,  No.  12,  pp.  2424-2431,  Dec.  1992. 

[24]  H.  Xin,  D.  E.  Oates,  A.  C.  Anderson,  R.  F.  Slattery,  G.  Dresselhaus,  and 
M.S.  Dresselhaus,  "Comparison  of  power  dependence  of  microwave 
surface  resistance  of  unpattemed  and  patterned  YBCO  thin  films,"  IEEE 
Trans.  Microwave  Theor.  Tech.,  48,  pp.  1221-1226,  2000. 

[25]  S.K.  Tolpygo  et  al.,  Deep  sub-micron  stud-via  technology  of 
superconductor  VFSI  circuits,”  Supercond.  Sci.  Technol.,  vol.  27,  no.  2, 
Feb.  2014,  Art.  ID.  025016. 

[26]  S.K.  Tolpygo  et  al.,  “Deep  sub-micron  stud-via  technology  for 
superconductor  VFSI  circuits,”  J.  Phys.  Conf.  Ser.,  vol.  507,  2014,  Art. 
ID  042043. 


